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Detection of Shocked Atomic Gas in the Kleinmann-Low Nebula 

M.W. Werner 1 , M.K. Crawford 2 , R. Genzel 2 . D.J. Hollenbach 1 , 

C.H. Townes 2 , and Dan M. Watson 3 

ABSTRACT 

We have mapped the 63/xm i P l —* 3 P 2 fine structure line emission of neutral atomic oxygen at the 
center of the Orion nebula with a resolution of 30”. There are three main emission peaks. One is asso- 
ciated with the region of strongest thermal radio continuum radiation close to the Trapezium cluster, 
and probably arises at the interface between the H1I region and the dense Orion molecular cloud. The 
other two [01] line emission peaks, associated with the Kleinmann-Low nebula, are similar in both dis- 
tribution and velocity to those of the 2/xm S(l) line of molecular hydrogen and of the high-velocity 
wings of rotational CO emission. The [01] line profiles in the KL region are resolved at a resolution of 
45 km s" 1 and have an intrinsic width of about 50 km s -1 ,. while the line profile in the vicinity of the 
Trapezium is unresolved. We show that the [01] emission from the KL nebula can be produced in the 
shocked gas associated with the mass outflows in this region and is an important coolant of the shocked 
gas. 

Subject Headings: infrared: sources - infrared: spectra - nebulae: Orion nebula. 

I. Introduction 

The Kleinmann-Low (KL) nebula in Orion is the site of a violent mass outflow associated with 
the formation of massive stars (c.f., the review by Genzel and Downes 1982). Infrared emission lines 
from highly excited states of FT, OH, and CO have been detected from this region (cf„ Gautier et al. 
1976; Storey et al. 1981; Beck et al. 1982); this emission is probably excited in shocks occurring where 

dense shells driven by the mass outflow strike the ambient molecular cloud material. In this paper, we 
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report mapping of the central 2.5' of the Orion-KL region with 30" resolution and show spectra taken at 
45 km s' 1 resolution which together reveal emission in the atomic oxygen i P ] -+ i P 2 ground state fine 
structure transition associated with the KL nebula. This transition is expected to be the dominant cool- 
ing transition for dense (10 8 n > 10 3 ) atomic interstellar gas for temperatures between 100 and 5000 
K. 


II. Observations and Results 

jl Observations 

The observations were carried out with the 91cm telescope of the Kuiper Airborne Observatory in 
1983 February and 1983 October, using the UC Berkeley tandem Fabry-Perot spectrometer described by 
Storey, Watson and Townes (1980). The beam had a FWHM of 30", with an effective area of a 38" 
diameter disk. In February 1983, the instrumental profile was Lorentzian with a FWHM of 170 km s" 1 . 
A chopper spacing of 4' in azimuth was used (approximately east-west). Spectral scans were obtained at 
approximately forty points (typical integration time 1 minute per point) on a grid spaced by 20" in both 
right ascension and declination. Each spectral scan included the line and the adjacent continuum, so 
that maps of line and continuum intensity were produced simultaneously with good relative positioning. 
Relative positions across the map are accurate to 5”. Absolute positions were obtained by assuming 
that the continuum peak is at R.A. = 5 / '32 m 46?8, Dec. = — 05°24'25" (1950, see Wynn-Williams et al. 
1984). Repeated measurements of the line intensity at a point near the map center agreed to ±10%, 
and it is this reproducibility, rather than the signal-to-noise ratio, which is the limiting uncertainty in 
the quality of the map presented in Figure 1. Observations were made during two flights in October, 
1983, at velocity resolutions of 60 and 45 km s -! towards several points on the map made in February; 
the 45 km s _1 spectra are shown in Figure 2. The beam size was again 30", and the chopper spacing 
was 5' in azimuth, which was again approximately east-west. Radial velocities were calibrated by 
measuring the relative wavelength separation between the [OI] lines in Orion (rest wavelength 
63.18372/tun) and an H 2 0 line at 63.3236/um. The line fluxes in Figures 1 and 2 were calibrated by 
measurements of the line to continuum ratio toward the continuum peak. We assumed a flux density 
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of 7 x 10 4 Jy for the KL continuum in a 30” beam (scaling by 0.6 the flux densities given by Werner et 
al. 1976). The overall calibration uncertainty is estimated to be ±25%. 

2 ) _ Spatial Distribution 

Figure 1 shows the spatial distribution of 63jum line and continuum flux over the region mapped, 
compared with other relevant data. The [01] line emission shows a peak 10” SW of the Trapezium 
cluster and a double-peaked ridge near the KL nebula. The peak near the Trapezium is very similar in 
position and shape to the 5 GHz radio continuum emission. [01] line emission in this region probably 
comes from dense (n M ~10 4-5 cm -3 ), warm (T ~a few hundred K) gas at the interface between the 
HII region and the molecular cloud. This interface region will be discussed in a forthcoming paper 
(Ellis and Werner, 1984; see also Jaffe and Pankonin, 1978). The [01] peaks in the KL region are 
close to but somewhat interior to the emission lobes of hot molecular hydrogen (Figure lb) and are 
also similar to the distribution of high-velocity rotational CO emission (cf. Erickson et al. 1982). The 
lower resolution data of Naylor et al. (1982) and Ellis and Werner (1984) also clearly show an exten- 
sion of the 63ju,m emission from the Trapezium toward the BN-KL region but lack sufficient resolution 
to identify a separate peak. 

3) Kinematics 

Figure 2 shows [01] spectra taken with a resolution of 45 km s -1 toward the radio peak west of 
the Trapezium, H 2 peaks 1 and 2 and the KL peak (the center of the 63/um continuum emission). The 
line toward the radio peak is unresolved; the line shape is identical with the instrumental profile. The 
LSR velocity of the line is 8.5 ±3 km s _l which is in good agreement with the velocity centroid of the 
6300A line of [01] at this position (Munch and Taylor, 1974). The line profiles near the KL nebula are 
resolved and the lines are significantly blueshifted. The line width at H 2 peak 1 is 67 ± 5 km s _l 
FWHM, and the peak of the line is at v LSR = 0±3 km s -1 . The line widths at H 2 peak 2 and the KL 
peak are 61 ±5 km s -1 and they are centered at v ISR = +3 ±3 km s _I . Most of the difference between 
these line profiles and the unresolved line toward the radio peak is probably due to the influence of 
high velocity [01] emission, blueshifted by about 20 to 50 km s' 1 . The presence of strong blueshifted 
wings is also characteristic of the IF 2/^m lines toward peaks 1 and 2 (Nadeau et al. 1982). The 
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intrinsic width of the [01] line towards the KL region may be estimated, by comparing the line width 
observed at 45 km s _1 with that at 60 km s" 1 resolution, to be about FWHM 50 km s" ! at Hi peak 1 
and 40 km s _1 toward the KL peak. These are lower limits to the width of the emission associated with 
the KL nebula itself since the profile probably includes a narrow component of line emission from the 
interface [01] gas in front of KL. 

III. Discussion 

_1) Properties of the [01]- Emitting Gas in the KL-nebula 

The [Oi]-emitting gas in the KL nebula has a large velocity dispersion and its emission has a spa- 
tial distribution similar to that of the 2/xm emission from hot hydrogen. It is, therefore, likely that the 
63/xm [01] emission is produced in shocks due to the mass outflow from the center of the KL nebula. 
The [01] line is of substantial importance to the gas cooling in the shocked region since its luminosity 
is of the same order of magnitude as the luminosities of the H 2 and CO lines. The peak surface bright- 
ness of the line near H 2 peak 1 is 3.3 ± 1 x 1CT 2 erg s~' cm" 2 ster~*, after subtraction of a background 
of about 2.5 x 1(T 2 erg s _1 cm" 2 ster -1 which may originate in the interface to the HII region. The 
total 63/xm line luminosity from the shocked region is 50 ± 25 L 0 , which is only about 3-4 times smaller 
than the total luminosity of all Hi infrared lines (—200 L 0 , c.f., Scoville et al. 1982, Beckwith et al. 

1983) , or the luminosity of the far-IR CO lines (150 L 0 ). The minimum column density of [01] gas 
required to produce the observed surface brightness is 4.5 x 10 17 cm -2 , assuming statistical equilibrium 
at high density (n^5 x lO 5 cm" 3 ) and T >> 230 K. This is probably a good assumption for the 
parameters of the [01] emitting gas in the KL nebula. The ratio of 63/xm intensity to that of the 
second [01] line at 145/xm toward the center of the KL nebula is between 20 and 40 (Ellis and Werner, 

1984) . This range is consistent with optically thin emission from dense gas (n H — 10'- 1 cm" 3 ) at high 
temperature (T —300 to 1500 K) and is also characteristic of shocked gas. With a width of 50 km s“\ 
the line has negligible opacity (r (l , — 10~ 2 ). Assuming all oxygen atoms are in the form of 01 and a 
cosmic abundance of oxygen (0/ II = 6.6 x 10“ 4 ), the lower limit to the column density of hydrogen 
nuclei in the [01] emitting region is 7 x 10 : ° cm" 2 , which corresponds to a mass of —0.15 M e . This is 
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significantly larger than the amount of very hot hydrogen gas behind the molecular shock (T >2000 K; 
Nji —5 x 10 19 cm -2 , c.f.. Beckwith et al. 1983) but less than the hydrogen column density estimated 
for the medium temperature molecular gas in this region (T —500— >1000 K; N|| —2 to 6 x 10 :| cm” 2 , 
cf. Watson et al. 1984). 

2) [01] Emission from the Cloud Shock 

The mass outflow in the KL nebula drives shocked, neutral gas into the ambient interstellar 
medium. As pointed out by Castor, McCray, and Weaver (1975), Hollenbach (1982), and Chernoff, 
Hollenbach, and McKee (1982, hereafter CHM), two separate shock fronts will exist under these cir- 
cumstances. The "wind shock" occurs at the point where the stellar wind collides with the inner edge of 
the shocked, expanding material. A second shock, (the "cloud shock") occurs where the outer edge of 
the expanding material strikes the ambient molecular cloud. Draine and Roberge (1982) and CHM 
have modeled the cloud shock as a C-type shock in which the magnetic field effectively "cushions" the 
shock front. These models have successfully reproduced the intensities of the high excitation molecu- 
lar lines of H 2 , OH and CO observed in the KL region. In such a shock, the [OI] 63/xm line is 
predicted to be two or more orders of magnitude fainter than the integrated H 2 emission, because any 
preshock atomic oxygen is rapidly converted to H 2 0 in the hot post-shock region. By contrast, the 
present observations show that the 63/um line is within an order of magnitude as intense as the total H 2 
emission. The column density of medium temperature (T >500 K) C-shocked gas, however, is 
sufficient to produce the [OI] emission if —10% of the oxygen could be maintained as [OI]. Thus, if 
our present understanding of the chemistry in this region is inadequate, it is possible that the [OI] 
emission could in fact be produced in the molecular shock. 

ll The Wind Shock 

Alternatively, the [OI] emission may come from the wind shock (Hollenbach. 1984; Draine and 
Roberge, 1982). If the outflow has a characteristic speed in excess of 100 km s' 1 , it strikes the inner 
surface of the compressed shell at a speed in excess of 50 km s' 1 . This is sufficient to dissociate and 
significantly ionize the wind material, which leads to formation of a J-type w-ind shock front with sharp 
density and temperature discontinuities and high post-shock temperatures (T »10 4 K). As the gas 
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behind the shock front cools and recombines, it passes through a range of temperatures from approxi- 
mately 5000 to 100 K in which [Ol] 63/irn emission is the dominant coolant for the largely atomic gas. 
Because the 63/a.m line is the dominant coolant in this temperature range, about 5000 k of energy (k = 
1.38 x 10~ IA erg K -1 ) is radiated in the line for each nucleon which crosses the shock front. Thus the 
flux of material into the shocks, n H v (cm~ 2 see -1 ’, is related to the surface brightness, I, of the 63/j.m 
line by: 

I = 5000 -p- n H v erg s _1 cm -2 ster -1 , (1). 

47T 

where n [( [cm -3 ] is the preshock hydrogen density in the wind, and v is the relative velocity between 
wind and shell. Detailed calculations over the range 40 km s~* < v < 200 km s _1 and 10 cm -3 < n H 
< 10 +;> cm -3 show that the [OI] 63/zm intensity lies within a factor of two of this linear relationship for 
n H v <10 12 cm -2 s~* (Hollenbach, 1984). Assuming that the wind shock in the Orion-KL region falls 
within this parameter range, the observed [Ol] 63^m intensity corresponds to a value of n H v = (6 ± 
2) x 10" cm -2 s -1 . In this limit, the total luminosity of the 63/um line gives a direct estimate of the 
mass loss rate M from the central object and is about 7 L 0 if M == 10 Mg y . For the observed total 
[OI] luminosity of 50± 25 L Q , the corresponding mass loss rate is (7 ± 3.5) x 10~ 3 Mg y _1 , which is 
comparable with or somewhat higher than other estimates of the mass loss rate in the KL region. If 
the [OI] emission arises predominantly in the wind shock, while the and CO infrared emission lines 
arise in the molecular shock, differences in the details of the line profiles and spatial distribution are 
expected. The [OI] data presented above are consistent with these expectations but do not exclude the 
possibility that the [OI] emission comes from the cloud shock. The [OI] 63/u.m emission line may be 
one of the few direct observational probes of J type wind shocks in embedded infrared sources. If [OI] 
emission from the wind shock can be unambiguously identified, its intensity may be an important meas- 
ure for determining mass loss rates in outflow sources. 
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Figure Captions 

Figure la (left): Contours of integrated 63/am [01] line emission toward the Orion-KL region with a 
resolution of 30" (heavy lines). Contour units are 2.5 x 10 -17 W cm -2 
(9.4 x 10~ 3 erg s -1 cm -2 ster -1 ). The dashed lines are contours of 5 GHz radio continuum emission at 
a resolution of 28" (from Johnston et al. 1983, contour unit 55 mJy per beam). The crosses mark the 
positions of the 63/u.m continuum peak (the KL-nebula, upper) and of 0'C Ori (lower), the brightest 
star in the Trapezium cluster. 

Figure lb (middle): Contours of integrated [01] emission (heavy lines), superposed on contours of v 
= 1 — * 0 (S ( 1 ) emission from H 2 at 2/zm (Gautier 1979), 30" beam, contour unit 
6 x 10~ 4 erg s -1 cm -2 ster -1 ). The northwestern, stronger peak is H 2 peak 1, the southeastern peak is 
H 2 peak 2. 

Figure lc (right): Contours of 63/zm continuum emission (heavy contours 1, 2, 4, 6, 8, 10, 12, 14, 16, 
18 times 3.8 x 10 3 Jy or 1.6 x 10~ 17 W cm' 2 ), superposed on the 5 GHz radio continuum map. 

Figure 2: Spectra of 63/am [01] emission toward several positions in the Orion-Nebula region: 10” west 
of Trapezium (bottom), H 2 peak 2 (second from bottom), 63/am continuum peak (third from bottom) 
and H 2 peak 1 (top). Beam size is 30" and the instrumental profile is a Lorentzian of FWHM 45 km 
s -! (indicated by thin lines in the upper three spectra). The bottom line is unresolved. The curvature 
in the continuum on the left and right of the [01] line is due to telluric H 2 0 absorption. The vertical 
scale is given separately on the y-axis for each spectrum. Typical integration times are 6 minutes per 
spectrum. 
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